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ABSTRACT 

A parameter-plane method, based upon a quasi-frozen 
assumption, for the design of electrical filters containing 
slowly varying nonlinear elements is presented together 
with examples of its application to the lowpass, highpass, 
bandpass and band-reject filter types. The nonlinear mag- 
nitude and frequency characteristics are presented in the 
form of parameter-plane graphs which allow value judgements 
to be made concerning the choice of nonlinear elements used 
in the filter design. The computer programs used to gener- 
ate the graphs are also included. 
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I. INTRODUCTION 



A filter is a frequency-sensitive unit which discrimin- 
ates between desired and undesired frequencies. Ideally 
the filter would have a unity-gain passband and a zero- 
gain stopband with no transition region between the bands. 
Also, its phase shift would be linear within the passband. 
Such filters do not exist today. The ideal filter is the 
optimum case and serves as a standard of comparison for 
all other filters. Filter design usually stresses either 
the magnitude response or the phase response depending on 
the intended use of the filter. 

Filters may be grouped into four distinct types based 
on their frequency discrimination: lowpass, highpass, 

bandpass, and band-reject. While classical design used 
passive, lumped, linear elements and image-parameter methods 
to realize these filter types, modern design uses the poly- 
nomial method of network synthesis and considers all types 
of elements in both their active and passive forms. The 
addition of nonlinear elements in these filters allows 
a closer approach to the ideal response than possible if 
the elements were linear. At the same time however, the 
design problem becomes more complicated. 

A parameter-plane method for the design of filter 
circuits is presented in this thesis. The method provides 
a relatively simple, easy-to-use way of incorporating 
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nonlinear elements into the design of electrical filter 
circuits . 
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II. THE PARAMETER PLANE METHOD 



The mathematical analysis of a linear dynamic system 
usually evolves into the evaluation of the roots of a poly- 
nomial formed from the differential equation description of 
the system. The root locations of this characteristic 
equation are then interpreted with respect to the system 
response characteristics. Any change in any circuit 
parameter value, however, changes the root locations. Some 
changes will improve the circuit response while others will 
degrade it . 

Mitrovic, Thaler, Siljak and others developed analysis 
methods which consider such parameter value variations [1] . 
The original method of Mitrovic was rather limited since 
only the last two coefficients of the characteristic 
equation were allowed to vary. Lines of constant damping 
ratio (c) were mapped from the s-plane into the plane of 
the variable coefficients using the characteristic equation 
as the mapping function. These transformed curves were 
then used to determine the variable coefficient values 
required to place the characteristic equation roots at 
desired locations. 

Generalization and extension of Mitrovic 's coefficient 
plane by Siljak and Hollister [2] led to the concept of the 
a, 3 parameter plane. This method considers the variable 
coefficients of the characteristic equation as nonlinear 
functions of two system parameters, a and 3. 
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F(s) 



n 

Z a. = 0 , a. = f . (a, 3) . 
i = 0 



Curve families of constant damping ratio and undamped nat- 
ural frequency are commonly mapped into the a, 3 parameter 
plane and used to infer system characteristics. 

In the case of a system transfer function, however, the 
characteristic equation is formed from the denominator. 

The numerator is not included in the formulation. There- 
fore, full frequency-response Information cannot be gained 
from using the characteristic equation in the parameter 
plane method. Thaler and Thompson developed a form of the 
transfer function which could be used in the parameter-plane 
method [3]. This transfer function form has been used by 
Glavls [4] and Staples [5] to make frequency-response plots 
and analyses of systems in which the coefficients contained 
linear combinations of a and 3. 

The derivation of the transfer function form developed 
in Ref. 3 follows: 

Let the system transfer function be given as 



_ N(s) 



where N(s) and D(s) are each of the form 



n 



N(s) = Z a.s^ 
1=0 ^ 
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and the terms are nonlinear function of two parameters, 
a and 3. The squared magnitude function can be obtained 
from 

m^ = T(jo))-T(-ja)). 

If the polynomials N(s) and D(s) are separated into their 
even and odd parts. 



T(s) 



N^(s) + N^(s) 
D^(s) + D^(s) ’ 



N^(s) - N^(s) 

= D^CsVV D°(s) 

N^(s) 2 - N (s)^ 
and T(s)-T(-s) = , 2 ^ 

D,(s)2 - D^(s)2 



From this it may be seen that both the numerator aind denom- 
inator are even polynomials. When s = jw both numerator 
and denominator are even in w. Thus 



Z 

S A . 0) 

rr,2 - i=0 ^ 

m = 



2i 



n 

S B, 0) 
k=0 



2k 



2 2 2 2 

where = f^ ( 1 , a, 3 , a3 , a ,3 ,ot 3,a3 ,...) 

and similarly for the terms. 

There are now four parameters present, alpha, beta, 
magnitude and omega. Any one may be expressed in terms of 
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the others and curves drawn depicting the various relation- 
ships among the parameters as defined by the transfer 
function. Phase versus omega may also be examined by 
rationalizing T(jco) and using the arctan relationship to 
solve for the angle. 

Although the original motivation for the parameter plane 
was to analyze linear dynamic systems, the method may also 
be used to analyze nonlinear systems. Consider a linear 
system. One of its characteristics is that a sinusoidal 
input signal of a given frequency will be repeated through- 
out the system as a sinusoidal signal of the same frequency. 
This is not always true in a circuit containing nonlinear 
elements. In most cases though, the internal signals are 
at least periodic and have the same period as the input 
signal. These signals may therefore be represented by a 
Fourier series. Thus a frequency-response transfer function 
may be defined relating the fundamental component of the 
output signal to the input signal and may be used in the 
parameter-plane analysis of a nonlinear system. The non- 
linearity is considered a parameter whose characteristic is 
defined by a locus of points in a parameter plane e.g., 
alpha versus transfer-function magnitude. The function 
a(v) is assumed to vary slowly enough that it may be con- 
sidered as "frozen" or fixed over the frequency range of 
interest . 

Further, this parameter plane analysis logic may be 
reversed and the method used as a tool in the design of a 
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system. Thus a desired system response may be specified 
and the a, 3 parameter variations necessary to produce this 
response determined by the method. 

The parameter-plane method would be far too laborious 
if done by hand. Therefore, FORTRAN IV computer programs 
have been developed by the author which allow the frequency- 
response form of the parameter plane method to be used in 
the design of filters with nonlinear elements. The program 
details are contained in Appendix A. 
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III. DESIGN OP FILTERS WITH NONLINEAR ELEMENTS 



Filters with nonlinear elements may be .designed by 
reference to parameter-plane plots of the voltage transfer 
function of the filter. The design may emphasize the mag- 
nitude response, the phase response or both. The nonlinear 
elements are represented as parameters alpha and/or beta 
and used in conjunction with the parameters of magnitude 
and of normalized radian frequency, omega. A computer 
generates selected plots of any one parameter versus another 
with a third parameter held constant and the fourth incre- 
mented to form a family of curves. A fifth parameter, 
phase, may be plotted versus omega for family parameters 
alpha and beta. The desired filter response is noted on 
these plots and the required characteristics of a particular 
nonlinear element determined. 

Prom given specifications a particular type of filter 
circuit is chosen and the transfer function determined. 

For simplicity let the nonlinear element be represented by 
the alpha parameter and be a function of output magnitude. 
Also assume that in this case the magnitude response is 
important while the phase response is not. 

The computer is told to make a Bode plot of the linear 
response of the filter for each of a family of alpha values 
and the response required by the specifications is over- 
plotted by hand onto the same graph (see Figure 1). Each 
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alpha-omega coordinate point along the specification curve 
indicates what the value of alpha must be in order to 
satisfy the specifications. These coordinate points are 
then picked off and transferred by hand to a computer- 
generated parameter-plane graph of alpha versus magnitude 
where beta has been held constant and omega incremented to 
form the curve family (see Figure 2). The locus of hand- 
plotted points defines the required characteristics of 
the particular nonlinear element with respect to the circuit 
output magnitude. If the nonlinear element were a function 
of frequency, the coordinate points would have been hand 
transferred to a computer-generated parameter-plane graph 
of alpha versus omega with magnitude as the family parameter. 
This would determine the necessary frequency characteristics 
of the nonlinear element. 

Filters with two or more nonlinear elements may require 
plotting points in the alpha-beta parameter plane (magnitude 
held constant and omega incremented or vice versa) to 
determine possible value trade-offs between the nonlinear 
elements. The curve plotted on the graph of parameter ver- 
sus magnitude or omega would still be used to define the 
nonlinear element characteristics. These same procedures 
could be followed with the response emphasis on phase rather 
than on magnitude. 

An additional consideration is necessary. The transfer 
function is the ratio of one voltage to another, E2/E1. 
Consequently, the plotted response is that of E2 with El 
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set to one volt. Thus in realizing the filter design the 
nonlinear element must be made a function of E2/E1 and not 
a function of E2 alone. A circuit to perform this division 
may be omitted if the nonlinear element is a function of 
frequency or when simulation sensitivity checks indicate 
that the filter response will be within acceptable limits 
for the predicted range of input values. 

Although the filters considered here are shown in a 
lumped R - L - C configuration, the R, L and C should be 
looked upon as representing voltage-current relationships 
rather than actual physical components. The lumped R-L-C 
could be transformed into distributed values and/or L - C 
interchanges made through the use of active networks such 
as gyrators, mutators, rotators, etc. Time invariance is 
assumed however. 

The steps in the design of a filter are: 

1. Form the transfer function, T(s), of the desired 
filter type designating the nonlinear elements as alpha, 
beta . 

2. Determine which computer-generated plots will be 
needed. 

3 . Encode data and use the computer programs described 
in the appendix to generate the plots. 

4. Select the Bode and/or phase plot as appropriate 
and draw in by hand the specified response curve. 
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5. Hand transfer the coordinate points along the res- 
ponse curve to the computer-generated parameter-versus- 
magnitude or parameter-versus-omega graph.s , 

6. Determine the characteristics of the nonlinear 
element from the locus of points plotted in step five above 

7. Check the design by computer simulation using 
several voltage input values to determine the effect of a 
changing input . 

8. If the nonlinear element is a function of magnitude 
decide whether or not the divider circuit is necessary. 

9. Frequency shift the component values from the 
normalized radian frequency to the frequency of interest. 
Impedance shifting may also be used. 

In the case of two or more nonlinear elements, step five 
will include the alpha- versus-beta plot for possible 
parameter value trade-offs. 

A. FILTERS WITH ONE NONLINEAR ELEMENT 

The following filter examples cover the four frequency 
types, lowpass, highpass , bandpass and band-reject. Each 
of these basic filters is made up of three elements, one R, 
one L and one C with the nonlinear element being a function 
of the output magnitude. The R element has been chosen as 
the nonlinear element since it may be readily associated 
with a field-effect transistor. The L or the C could just 
as well have been chosen. 
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1 . Lowpass Filter 



Figure 3 illustrates one possible lowpass filter 
configuration. Since frequency scaling will be employed as 
the last design step, the L and C values are normalized to 
one. The parameter alpha replaces the nonlinear element R. 
Thus the circuit transfer function is 

T(s) = i y . 

1 + as + s'^ 

The necessary computer plots for the design of this 
circuit are the magnitude-versus-omega (Bode) and alpha- 
versus-magnitude plots. 




Figure 3. Lowpass Filter 



Figure 4 is the Bode plot for the circuit. The alpha 
parameter curves on the plot indicate the output magnitude 
of the filter for the normalized radian frequency range of 
0.1 to 10.0 radians when alpha is fixed at a particular 
parameter value. Thus if alpha were set equal to 1.0, the 
output-versus-frequency response of the filter would be 
that indicated by the curve for alpha equal to 1.0. There- 
fore, this graph illustrates the response of a linear filter 



21 




NICHOLS. BGDf PLOT, PPRRMErE-;R--fiLPHq 
LONPRSS filter, L-C---1.0, R---RLPHR 



22 



OHEGR 



for several values of one of its elements, R. The filter 
specifications are indicated by the heavy line drawn on 
the graph by hand. In this case it is assumed to be the 
magnitude response of an ideal filter. The zero gain in 
the stopband of the ideal filter is represented by a line 
at the minus 60.0-dB level. 

The coordinate points of the intersection of the alpha 
parameter curves and the desired response line in Figure ^ 
are picked off and transferred by hand to the alpha- versus- 
magnitude parameter-plane plot in step five of the design 
procedure. This plot is shown in Figure 5- The heavy 
line is the locus of points transferred. It indicates the 
desired characteristics of the nonlinear element when the 
value of the element is a function of the filter output 
magnitude . 

The initial value of alpha for this plot is less than 
or equal to 2.0 at an omega of 0.1. The required value 
of alpha at cutoff is 1.0. For a flat, zero-dB response 
between these two omega points, required alpha values are 
1.5 at an omega of 0.2 and 1.35 at an omega equal to 0.5. 
The plot also indicates that an alpha value of 1.0 would 
deviate from the ideal response by +1.0 dB at omega equal 
to 0.5 and less than about +1.2 dB between 0.5 and cutoff. 
If this is tolerable, alpha may be fixed at 1.0 for values 
of omega less than 1.0 and the design problem simplified. 

In the region where omega is equal to 1.0 the alpha- 
magnitude relationship is logarithmic. The filter output 
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is 0.0 dB at alpha equal to 1.0, -60.0 dB at alpha equal 
to 1000 and - 80.0 dB at an alpha of 10,000. The minus 
60.0-dB level requires alpha to have a three-decade range 
in value. The minus 80.0-dB level is attainable with the 
four-decade range of field-effect transistors. Alpha may 
remain at this last value for frequencies above 1.0. If 
this is so, the output level will fall off at a 20-dB/ 
decade rate. 

At this point consideration should be given to the 
realization of the circuit. Although theoretically any 
function may be generated, the economics of the situation 
usually limits the degree of circuit complexity allowed 
and causes simplifying approximations to be investigated. 

The fixing of alpha at 1.0 for the frequency region below 
1.0 and to a fixed value for the region above 1.0 would be 
simplifications in the case of the lowpass filter. The 
drop in output magnitude below 0.0 dB at cutoff when alpha 
is equal to 1.0 may be used to initiate the required 
parameter value change. The change in R could come from 
the filter transfer-function output voltage being subtracted 
from a fixed reference level and the difference applied to 
the gate of a field-effect transistor. 

Figure 6 is the phase-plot companion to the Bode plot 
of Figure 4. The dark line indicates the expected phase 
response of the filter when alpha equals 1.0 at values of 
omega below 1.0 and equals 10,000 at omega values above 1.0. 
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FIGURE 6. 
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The phase fall-off in the latter region is very slov;. 
Designing for the ideal magnitude response does not yield 
a linear phase filter. 

If alpha were a function of frequency, the plot of 
alpha versus omega. Figure 7, would replace the alpha- 
versus-magnitude plot. Again, the heavy line is the locus 
of points transferred from the Bode plot of Figure 4. The 
0.0-dB curve defines the exact values of alpha necessary 
for a completely flat magnitude response in the passband 
region. The area between the 0.0 and +1.0 dB curves indi- 
cates the latitude in alpha values available for frequencies 
below 0.2 radians. 

The computer simulation and input sensitivity check 
results are depicted in Figure 8 for the output magnitude 
and in Figure 9 for the phase. A 1.0-dB region of uncer- 
tainty or lag is assumed for the output level at which the 
parameter value change is initiated. Alpha is fixed at 
1.0 for output values above the minus 1.0-dB level and 
allowed to go as far as 10,000 whenever the output drops 
below -1.0 dB. The curves represent the circuit response 
when the divider circuit is omitted and the input voltage 
varies around 1.0. 

2 . Hlghpass Filter 

Figure 10 illustrates a possible highpass filter 
configuration. As in the lowpass filter circuit L and C 
are normalized to one and alpha represents the nonlinear R 
in the transfer function 
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T(s) = 



1 + as + 



s 



2 • 



The necessary computer plots for the design of this filter 
are the Bode and the alpha- versus-magnitude plots 




Figure 10. Highpass Filter 




Figure 11 is the Bode plot for this circuit. Since the 
circuit configuration is related to that of the lowpass 
filter by an interchange of frequency-sensitive elements, 
the Bode plot is a reflection of the lowpass filter circuit 
about the omega equal to 1.0 axis. This relationship is 
further brought out by the alpha- versus-magnitude plot of 
Figure 12. It is exactly the same as Figure 5 except for 
reversed omega values. 

Circuit realization is also related to that of the low- 
pass filter. In this case the parameter alpha could be 
fixed at 10,000 in the omega region below 1.0 and allowed 
to decrease to 1.0 in the region above 1.0. The filter 
transfer-function output would be added to a reference level 
and the sum applied to the gate of a field-effect transistor. 
The change in alpha would be initiated when the output mag- 
nitude exceeded a given value, say -80.0 dB. 
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3 . Bandpass Filter 



The three-element circuit configuration for a band- 
pass filter is shown in Figure 13. Again L and C are set 
to one and R replaced by alpha. The transfer function is 



1 + as + s^ 

The Bode plot and the alpha-versus-magnitude plot are the 
computer plots used in the filter design. 




The heavy line drawn on the Bode plot of Figure l4 
specifies a 1.0-omega-wide bandpass symmetrical about the 
point where omega equals 1.0 and a stopband level of 
-60,0 dB. 

Figure 15 shows the specifications in terms of the 
nonlinear element versus magnitude. The desired alpha 
characteristics in this case are more complex than in the 
lowpass case. A small change in alpha causes a large 
output change when alpha, itself, is small. When alpha is 
large, which is the case in the passband, large parmeter 
value changes are necessary to affect the output. This 
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would indicate a relatively unstable attenuation level in 
the stopband region but a stable level in the passband. 

The width of the passband is limited by the minimum 
attenuation required in the stopband. The wider the pass- 
band, the higher the maximum alpha value must be for a flat 
response. The range of values available with the chosen 
element would then Indicate the point where the stopband 
attenuation is too low to meet specifications and thus 
limit the width of the passband. On the other hand, the 
filter passband could be made very narrow. 

Realization of the filter can not be done on the basis 
of magnitude alone. Some sort of frequency or phase 
sensing device is needed to determine the region of opera- 
tion, either above or below omega equal to 1.0. Given a 
frequency sensor the circuit could be formed by letting the 
magnitude response level crossing initiate the alpha 
parameter value changes to define the end points of the 
passband. This is the same type of operation as used in 
the lowpass and highpass filters. 

The phase plot is shown in Figure 16. The zero phase- 
shift region indicates the passband where alpha values are 
large. The phase below the passband would be 90 degrees 
leading and above the passband 90 degrees lagging. 

4 . Band-Reject Filter 

A three-element band-reject filter is illustrated 
in Figure 17. The transfer function when L and C are set 
to one and R replaced by alpha is 
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T(s) 



1 + s + 

The design of the filter calls for Bode and alpha-versus- 
magnitude plots. 




Figure 17* Band-Reject Filter 

The specifications are indicated on the Bode plot. 

Figure 18 , and the alpha-versus-magnitude plot. Figure 
19* The phase plot is shown in Figure 20. The plots 
confirm the expected, that the band-reject filter is the 
inverse of the bandpass filter. 

As in the case of the bandpass filter, the band- 
reject filter can not be realized by using magnitude res- 
ponse alone. It too must have a frequency or phase sensor 
to indicate the region of operation. 

B. FILTERS WITH MORE THAN ONE NONLINEAR ELEMENT 

More complex filters may be formed using more than one 
nonlinear element. They may be built up by combining filter 
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R.R. NICHOLS, BODE PLOT, PRRflHETER=fiLPHfi 
BRND-REJECT FILTER, L=C=1.0, R-RLPHR. 
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FIGURE 19. 

R.R-,HJCH0LE>.RLPHR U5 HRCNITLIDEIDB) .PRRRH-QMEGR 
BAND -REJECT FILTER, L--C---1.0. R=RLPHfl. 
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FIGURE 20. 

R.fi. NICHOLS. PHRSE PLOT. PRRRMETER=RLPHR 
BRND-REJECT FILTER, L-C--1.0. R=RLPHR. 
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sections and/or by using several nonlinear elements in the 
same section. 

1 . Concatenating Filter Sections 

Additional attenuation of undesired frequencies may 
be obtained by concatenating filter sections. The sections 
are designed as individual units and then coupled together 
by buffer amplifiers. 

A considerable improvement over the single-section 
design of a bandpass filter may be made by concatenating 
lowpass and highpass sections. The passband end points 
are set in the frequency shifting step by the amount of 
HP-LP cutoff frequency overlapping used. A band-reject 
filter may be designed by inverting the bandpass filter 
output and adding it to the input signal. 

2 . Cascading Filter Sections 

There are two types of cascaded filters. The first 
type consists of those in which the nonlinear element in 
each section is made a function of the transfer ratio of 
the total filter. This type also includes those filters in 
which the nonlinearity is a function of frequency since 
this parameter, or at least its first Fourier approximation, 
is common to the filter as a whole. 

Figure 21 shows an example of this first type of 
cascaded filter. Two lowpass filter sections are cascaded 
with both R1 and R2 functions of the total filter output 
magnitude. R1 equals k times R2 . Setting L, C and k to 



^3 



one and replacing the R’s with alpha gives a transfer func- 
tion of 



2as + (2 + a^)s^ + 2as^ + s 




This type of cascaded filter is designed as a single 
filter since the transfer function contains only one param- 
eter, alpha and it is a function of magnitude. The Bode 
plot for this circuit is shown in Figure 22, the alpha- 
versus-magnitude plot in Figure 23 and the phase plot in 
Figure 24. The specifications are assumed to be the same 
as those for the single-sections lowpass filter except 
that the stopband attenuation is now set to -120.0 dB. 

In this filter design there is no simple approximation 
which would give a reasonably flat response below cutoff. 

The alpha value must change smoothly from 5 • 0 at omega equal 
to 0.1 to- 0.1 at omega equal to 1.0, The relationship of 
alpha to magnitude in the omega equal to 1,0 cutoff region 
is shown in Figure 23. This is not a simple logarithmic 
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R. ft. NICHOLS. BODE PLOT. PflRflhETER=flLPHR-R 
2 -lP sections CflSCflDEO. TVPE 1. L---C-1»0 
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R.q.NJCHOLS-nLPHfl(R) US NflCf DB) , PflRRM=aMEGfl 
2-LP sections CflSCRDED, TYPE 1, L=C--1.0 
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FIGURE 24. 

R.R. NICHOLS,. PHRSE PLOT. PflRRMETER=flLPHfll R) 
2-LP SECTIONS CflSCRDED. tvpe 1. L--C--1.0 
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relationship as was the case with the single-section low- 
pass filter. 

This type-1 cascaded, lowpass filter has an advantage 
over the two-section concatenated filter in that it does 
not need the buffer amplifier; but, on the other hand, it 
does require a more complex parameter characteristic in 
the passband. 

The second type of cascaded filter is illustrated in 
Figure 25. It consists of those filters formed by the 
direct addition of sections. In forming the transfer 
function L and C are set to one, R1 is replaced by alpha 
and R2 by beta. The transfer function for the total filter 
is 



T(S) = E3/E1 = p p jr . 

(a+8)s + (2+a8)s + (a+8)s'^ + s 

R2 is a function of E3/E2. The transfer function for this 
second section is 

» 

T(s) = E3/E2 = p . 

1 + 8s + s 
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In addition, R1 is a function of E2/E1 so that 



T(s) = E2/E1 = -r s ^ 

(a+6)s + (2+a3)s'^ + (a+6)s'^ + s 

A Bode plot for the E3/E1 transfer function is shown 
in Figure 26. The parameter is alpha with beta set to 1.0. 
The specification curve is drawn in as in the previous 
examples with the stopband gain level at -120.0 dB. Due to 
the alpha/beta symmetry in the transfer function. Figure 26 
is also valid for beta as the parameter with alpha set to 
1.0. 

R1 is a function of an internal voltage-transfer ratio, 
E2/E1. The method of design in this case may be illustrated 
by assuming that R2 is a fixed value and that R1 is the only 
nonlinear element in the circuit. As before, the alpha-omega 
coordinate points of the desired response curve are picked 
off from the Bode plot and transferred to the alpha-versus- 
magnitude plot illustrated in Figure 27 and the required R1 
characteristics determined. The only difference here is 
that the magnitude used is the E2/E1 voltage ratio rather 
than the output magnitude of the E3/E1 transfer function. 

In this filter, however, R2 is also nonlinear so the 
plot of R2 versus E3/E2 magnitude must be considered. This 
plot is the same as that of the single-section lowpass 
filter in Figure 5* Therefore, realization of the 'passband 
region may be simplified by setting R2 at 1.0 and causing R1 
to follow the characteristic curve in Figure 27. The values 
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R.R. NICHOLS. eODE PLOT, PflRflHETER=flLPHfi=Rl 
(E3^EI).2-LP sections CRSCHDED. R2--BETfl-l, 0=L-C. 



50 



OMEGR 




R. n.lilUiJl h'LRHHlH" 
! b 4 ’ t ■ ' = 4' I JLt, I 



Mhu' IJtJ) .RHRhM^JMLRH 
LSbUHDtLi.h?- 1 ^ L- ’ . 0 



■51 



rpc 



of alpha to be used in the passband region are determined 
from Figure 28. This figure is a plot of alpha versus 
beta for passband frequencies where the filter output 
magnitude is 0.0 dB. 

In the cutoff region R1 and R2 change from their pass- 
band low values to their stopband high values. R2 
follows the omega equal to 1.0 logarithmic curve to the minus 
80.0-dB level at beta equal to 10,000. A four-decade range 
R1 element goes from 0.1 at 0.0 dB to 1000 at -60.0 dB . In 
so doing it follows the hand-drawn curve in Figure 29. This 
curve is a composite of the line for omega equal to 1.0 in 
Figure 27 where R2 equals 1.0 and that of Figure 29 where 
R2 equals 10,000, The chosen values are determined from 
the E2/E1 alpha- versus-beta plot of Figure 30 where omega 
is fixed at 1.0 and output magnitude is the curve family 
parameter. 

In the stopband region R1 may be left at 1000 and R2 
at 10,000. The resultant attenuation at cutoff is -1^0.0 
dB from which it falls off at ^0 dB/decade. Figure 31 is 
the Bode plot of E3/E1 for the stopband region above omega 
equal to 1.0. 

Figure 32 illustrates the phase response for passband 
values of alpha and beta while Figure 33 shows the phase 
response for stopband values. 

Figure 3^ and 35 show what the output magnitude and 
phase response would have been if two lowpass circuits 
had been added together directly. The response is not the 
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FIGURE 28. 

R.R. NICHOLS. RLPHfl(Rl) US BETfi(R2) .PRRRM=OHEGfl 
2-LP SECT CflSCflDED.L-.^C=1.0.hflG(DB)^0.0 
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L^fl'^HICHOLS.flLPHR(Rl) US 

(t2/ElU2-LP SECT Cfl5CRDED.L-=C--1..0,0HEGfl -UO- 
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R.H.HJLHdLS, BODE PLOT, PHRRHET[ R=BE I 
iti.El.U2'LP SECT CRSCHDED.Rl-lOnO..O»l,-l>-,1.n 
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FIGURE 32. 

R.R. NICHOLS. PHRSE PLOT. PRRRHETER-BETR=R2 
(E3/E13.2-LP SECT CR5CRDED.R1=RLPHR=1. 0-L=C 
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R.R. NICHOLS. PHRSf PLOT. PRRRHETLR-=BETq=R2 
I ES-El) , 2-LP SECT CnSCRDED.Rl=1000.G.L=-C-1.0 
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FIGWRE 34. 

R.fl.NlCHOLS.BQDE Pl0T.PRRRM-(UINPUT*1.0) 
2-LP SECT CRSCnOED.TVPE 2.SIMULRTiGN CHECK 
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R . n . N ] CHOLS , PHASE PLOT , PARRH= ( U I NPUT*1 . 0 ) 
2-LP SECT CASCADED. TYPE 2. SIMULATION CHECK 
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same as would have been obtained by concatenating the two 
sections. The concatenated filter is easier to realize 
than a type-two cascaded filter and has a greater stopband 
attenuation . 

3 • Two Nonlinear Elements in the Same Section 

A third configuration for two nonlinear elements 
may be obtained by letting both the R and the C in the basic 
lowpass filter be nonlinear elements and functions of the 
filter transfer ratio. The design procedure is the same 
as that of the previous sections. Bode plots are made for 
each of the nonlinear elements with the other set to a best 
guess value. Prom the specifications and the Bode plot 
information minimum and maximum values for the elements are 
chosen to use in a design trial. The information is 
transferred to parameter- versus-magnitude plots and the 
nonlinear element characteristics determined. The inter- 
dependence of R and C requires a composite curve to be used 
in defining these characteristics. 

Figure 36 is the alpha-versus-bet a plot for omega equal 
to 1.0. Alpha represents the nonlinear R which is assumed 
to have a four-decade range of values. Beta represents the 
nonlinear C which is assumed to have a two-decade value 
range. The specifications are those of a lowpass filter. 

The heavy line in Figure 36 indicates the chosen alpha/beta 
relationship in the filter cutoff region. Figures 37 and 
38 are the alpha- and beta- versus-magnitude plots with the 
element characteristic curves drawn in. Computer simulation 
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pfCURF 36. 

R.q.HICHaLS.qLPhR 05 SETR.pRRRHETER-HRC'DSI 
lP f [lTER.. l-v1. 0. OHEGR:--!. G, R--RlPHR, C -BE tR 
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FIGURE 3?. 

R.R.NICHOLS, ALPHA US HRCf DBUPflRRH=aNEGfl 
LOHPRSS FILTER, L-1. 0. R--RLPHR, C-BETfl=l. 0 
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FIGURE 38- 

R-fl- NICHOLS. BETH US NRG(DB) .PRRRHETER-OHEGR 
LOWPflSS FILTER.L--l-O.R-RLPHR=l-O.C--BETR 
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results for this design are shown in Figure 39 for output 
magnitude and in Figure 40 for the phase. Alpha and beta 
are set to 1.0 when the output magnitude is greater than 
- 1.0 dB and allowed to go to 10,000 and 100 , respectively, 
when the magnitude drops below - 1.0 dB . As before, the 
- 1.0 dB level is used to represent a lag in changeover for 
the change initiation at 0.0 dB. The cutoff attenuation 
is 120.0 dB while the phase response is nonlinear. 

C. DESIGNING FOR PHASE RESPONSE 

When the filter specifications concern phase response 
rather than output magnitude the desired response is 
drawn in on the phase plot rather than on the magnitude- 
versus-omega plot. The alpha-omega points are still 
transferred to the alpha-versus-magnitude plot for defini- 
tion of the characteristics of the nonlinear element when 
it is a function of the filter output transfer function. 
Where the magnitude of the output is of interest the same 
alpha-omega points may also be transferred to the magnltude- 
versus-omega plot for an indication of the magnitude response 
corresponding to the filter phase specifications. 

Figure 4l is the phase plot for the basic lowpass 
filter of Figure 3« The heavy line represents a specifica- 
tion calling for a linear l80-degree phase response between 
omega equal to 0.2 and 1.8. The alpha-versus-magnitude 
plot is shown in Figure 42. The magnitude characteristic 
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l^’lGURf; 40. 

R. R. NICHOLS = PHASE PLOT. PRRRM-( UINPUT*!. 01 
LP FILTER SIHULRTION CHECK. L--1.0.R--RLPHR.C--BETR 

. 67 







figure 41n 

R.fl, NICHOLS. PHflSf US OHEGR. PRRflH--RLPHfi-R 
PHRSE DESIGN HITH LOHPRSS FILTER. L-C-1-0 
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R.R.N]CHQLS,nLPHRlRJ US HflG( DB) , PRRRM-OhEGR 
PHASE DESIGN R]TH LOhPRSS FILTER, L--C=i. 0 
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alpha is fairly flat for the lower half of the frequency 
range but falls off in the upper half. This may also be 
seen in the magnitude-versus-omega plot of Figure 43. 

The basic lowpass filter may be designed for a 180- 
degree linear phase shift if the output magnitude is not 
part of the specifications. If, however, a flat magnitude 
response is also important, the maximum phase shift would 
be limited to 90 degrees for this frequency range. 
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R-R-NICHOLS.HRGfDB) US OHEGfl.PflRRH=flLPHfl--R 
• PHRSE DESIGN WITH LONPRSS FILTER. L-C-l-O 
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IV. CONCLUSIONS 



Filters with nonlinear elements may be. designed using 
the frequency response form of the parameter plane method. 
The method is a graphical one which might be thought of as 
a counterpart to the describing-function method used in 
analytical design. Parameter-plane plots furnish the mag- 
nitude and frequency characteristics of the nonlinear 
element so that the filter may be designed to meet specifi- 
cations. This method has the advantage of presenting the 
designer with a picture of the possible circuit conditions 
and thereby allowing value judgements to be made concerning 
the particular response chosen. 

On the other hand, this form of filter design has the 
disadvantage common to all graphical methods. As the num- 
ber of variables increases so does the number of graphs 
needed and the rate of growth of numbers of graphs is 
much faster than that for the number of variables. The 
point is finally reached where the number of graphs needed 
surpasses the ability of the designer to interpret and 
integrate them. At this point a hybrid form of graphical- 
analytics must be used. 

The basic building block for filters with nonlinear 
elements is the lowpass filter. With a nonlinear R having 
a four-decade range of values, it will provide a stopband 
attenuation of 89 dB and a passband with a maximum rise of 
less than 1.5 dB. Adding a nonlinear L or C with a 
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two-decade range of values will increase the stopband 
attenuation to 120 dB. Multiple filter sections may be 
concatenated to reach even higher levels of attenuation. 

If the L and C elements are interchanged, the lowpass 
filter is transformed into a highpass filter. A bandpass 
filter may be formed by concatenating low- and highpass 
sections. The band-reject filter may be generated from the 
bandpass form. The lowpass filter may also be designed 
for a linear phase response. A linear l80-degree phase 
response is available for cases where a flat output mag- 
nitude is not necessary, or a linear 90-degree response 
where a reasonably flat output magnitude is specified. 

Several areas exist for future investigation. The 
first one is that of the physical realization of these non- 
linear filters, primarily the lowpass filter. Another is 
the investigation into the response of the filter to 
complex wave shapes. A third area is the extension of the 
method to include more complicated circuit configurations 
where the question of the optimum mix of nonlinear elements 
does not have an apparent answer. In addition, the use of 
feedback and feedforward paths should be investigated as 
response-shaping methods. 
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APPENDIX A. COMPUTER PROGRAM INFORMATION 



The two computer programs are written In the common 
FORTRAN IV language for use on the IBM-36O model 67 computer 
with a CALCOMP model 765 plotter for graphics. This is the 
installation at the Naval Postgraduate School. Documenta- 
tion including data-card format information and internal 
section headings is written into each program. Both programs 
are capable of expansion beyond the present limitations of 
tenth-order transfer function and third-order alpha/beta 
parameters with I6 curves per graph and 900 points per 
curve. Linear, semi-log and log-log plots are available. 

I. PARTFI PROGRAM 

PARTFI provides graphs of alpha, beta, and omega (ordi- 
nate) versus magnitude (abscissa). Numerator and denomi- 
nator polynomials in omega squared are formed from the 
s-domain transfer function and magnitude squared determined 
from their quotient. Graph points are generated by incre- 
menting the desired ordinate parameter through the selected 
range of values while the other two parameters are held 
constant. The curve-family parameter is then set at the 
value for the next curve and the process repeated. Phase 
values are determined along with magnitude whenever the 
incremented parameter is omega. These phase values are 
available for plotting along the abscissa against omega. 
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Subroutines used with PARTFl are GRPLOT, ABVAR and 
INTTFN . GRPLOT provides the plotting package for the 
program. The other two are used with the circuit simula- 
tion feature. 

ABVAR is used to enter the FORTRAN statements describ- 
ing the nonlinear parameter characteristics. Alpha and 
beta are initialized in the main program at each incremen- 
tation of omega, a solution made for magnitude and the 
information passed to this subroutine. The values of 
alpha, beta, magnitude and omega are then available for 
logic testing to determine the required alpha and beta 
values for that circuit condition. Flag KNR will cause a 
recycling of the magnitude solution using the new values 
of alpha and beta until the change in alpha is less than 
EPSA and in beta less than EPSB. KKl is available as a 
loop counter. It is initialized whenever alpha and beta 
are initialized. 

In cases where one of the parameters is not a function 
of the output transfer function, INTTFN is used. It is 
called by ABVAR when MDATA is set to one in the data set. 

The subroutine reads in additional data describing the 
internal transfer function and uses the ABVAR information 
to compute magnitude. This value, X2, is then passed back 
to ABVAR for use in determining parameter values. 

The storage requirements of PARTFl are set primarily by 
the "DIMENSION PHASE (16,901), YY(l6,901)" card. The number 
"l6" represents the maximum number of curve parameters or 
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voltage input values used while "901" is the maximum number 
of computation points plus one. Using these values, the 
program requires 205K plus 9K for the plotting package 
and about 19K for operations. This comes to about 233K. 

This may be reduced considerably by changing the statement 
to reflect the actual number required. For instance, if 
only six curves were to be plotted, the storage requirements 
would be reduced by about 72K (2 X 10 X 901 X 4). If only 
500 points were used for each of the six curves, a further 
reduction in storage of about 19K would be gained. 

As expected, the number of curves and plotting points 
also affects the run time. Compile time is about 23 
seconds and the time for the link step about 1.5 seconds. 

A rough rule of thumb for go time is four seconds per 
900-polnt curve. A six-curve run of 900 points per curve 
vould have an execute time of around 48.5 seconds or about 
50 seconds if a phase plot were also made. 

Figure 44 is a repeat of Figure 6 using only 300 points 
per curve. Figure 45 uses only 100 points per curve. For 
design graphs of this size it would be wasteful in both 
storage and run time if 900 points per curve were demanded. 

2 . PARTF2 PROGRAM 

• PARTF2. provides graphs of alpha and beta versus omega, 
and alpha versus beta. As in the first program the input 
s-domain transfer function is separated into numerator and 
denominator functions which are individually transformed 
into the omega-squared domain. In PARTF2, however, the 
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FIGURE 44. 

R.R. NICHOLS. PHRSE PLOT, PRRRMETER-RLPHR 
LOHPRSS filter, L-C--1.0, NRPNTS--300 
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FIGURE 45. 

R.R. NICHOLS. PHASE PLOT. prrrhETER^RlPHR 
L0WPRS5 FILTER, L-C'-i-O. NRPNTS-iOO 
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denominator is then multiplied by the square of the magni- 
tude and subtracted from the numerator to form a single 
polynomial in terms of the unknown ordinate parameter. The 
real roots of this polynomial together with the incremented 
abscissa values generate the graph points for the curve. 

The subroutines used with PARTF2 are GRPLOT, P0LRT2 and 
POLRT 3 . GRPLOT provides the plotting package for the 
program. The other two are root-finding subroutines. 

P0LRT2 is used when the polynomial order is two and POLRT 3 
when the order is three or more. Exponential overflow and 
underflow may occur in the high-order terms of the poly- 
nomial in POLRT 3 . This will not affect accuracy. The 
ERRSET command cancels the overflow and underflow error 
print order and prevents the program from being prematurely 
ended when a given number of them have occurred. P0LRT3 is 
a modified POLRT routine from the IBM System /360 Scientific 
Subroutine Package. 

Contrary to PARTFl there is no one predominant factor 
setting the storage requirements. The limiting case of each 
factor contributes equally. Including the plotting package 
and operating space requirements, the storage requirement 
comes to about 148K. 

PARTF2 run times are much longer than those for PARTFl. 
Forming the polynomial, solving for the roots, selecting 
roots, etc. multiplies the run time over that for the simple 
quotient procedure of PARTFl by a factor of seven. The go 
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time rule of thumb is 28 seconds per 900-point curve, 
compile, link and edit time is about 23 seconds. Thus 
estimated run time for eight 900-point curves would be 
about 2Hj seconds (8 X 28 + 23). It pays to use fewer 
curves and fewer points with PARTP2 . 
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*DX,nY,7PMF , ALPHA Tc 



T'^265 



CALL ERPSt^T ( 207 , 300 , 300 , 2 , 1 ,? 0 «? ) tf 2662 

TFIT = 0 TC 7 A(t 3 



OOO 



xi u s^r^ ecu Or-<(NjO xj UA'C'r^o^u c^r-«fN r» ir^>or' a u o»-^rv^r^t>^u >ur* a/f»^Or-^(xr '-4 u 
^<)^0'X)^x. >O'0r"r' r^r'-r‘r^r r^r^r^cca’Q oca o^a cua. ecu u u u uu u u u u'ooooooo 

'41 >0'rvrND'X’'fWw c»Ni’>/’sxw 'J nT'O'U'^ 'r^;r r^r r^ 

fvjfMrvjfMc\j fNirv fMc\jrMrvjcvjCM(Mrvrv»fMrN‘rN)<\j(Mfvjcvj(N (njcvkmc'O' fvjr fNjrvjr fMf\ rvtrvjoiO'Tvjrvjfv 

ULiLLLCULiUUQU UUULtLi U UULiUUULU ULLLLULiUUULLUUULLL LiUaULiLL 







¥ 






¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥r 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 


¥ 




¥ 






¥ 


or 




¥ 


ijj 




¥ 


K 




¥ 


Z 




¥ 






< 


6 






o 




(/' 






IX 


o 




o 


z 




c? 


<r 




o 








i/' 




or 


UJ 




u 




o 


QL 


-j 




CL 


<3 




lU 


> 


in 






rvj 


K 


OK 




IL 


f^a. 


o 


ir 




fHCvl 


IT 




ro^ 


r'-fr 









H 

¥ 

it 

¥ 

-H 

¥ 

¥■ 

¥ 

¥ 

¥ 

> 



¥ 

¥ 

¥ 

¥ 

¥ 

< 

¥ 

r 

¥ 

¥ 

¥ 

¥ 

¥ 

¥ 









> 






xt 








t/' 


c 








z 






o 


d 






K 


o, 

ir^o X 






U 


•^O 






CD 










ir ►mu 




CC 


«-« 


►m ►u^ 




zz 


r- 


Oir.-^ 


r-H ^ 


►— »— t 


ai r-‘ ^ 




“5 UL 


^ w 


ooo 


o ►-{; 




LJL.1 


i • 






q:*-' 


> X .-*00 




4* X 


.-.XX 


¥ ¥ •OO 


.-•o *>»• 


<f II ^ II 


>0 -4 


OOU-U OU‘ 


4- U -X— '-t 


z-l 1 ^ 


xtuu aj 


• •■^C3^^xJ* 


1- 1 


I-4K 


Z-i-J 


OCJ *OCJ 


Cl-^ 1 ♦ 


II o-?> •- 


C^cLcr 


• •C' 


►-a;a z ♦ 



in¥ 

m¥ 

^ ¥ 

'-'U'.-k CJ O 

LLr-^^l f\i CM C\J^¥ r 

(D ¥ ^ r^'¥ 

^ycx.<i c mo<vo I ¥ : 

II II M K Ml- II K ^¥ 2 II II K- II QOOO I I XOOK- II CO*- 

i Q.I-Q: a CC w^llXIlrH IILLC CJ II II II X II II W II II K'-^'-HlX>ll 

II UJlCUL^f ULCUjOL OUiCU * XXC0->C»- C^O-^OC Z H CCU OOOOU U-LL ♦ II HZ 

I— cr*— o^-o^tr^-'ii zzzxos o*-*o*-»x>-i—xx>*—xxc*-»*-^ x>-^- 



-ft 




* 








¥ 


o-ff- U>CJU' 


O 


04 U' 


CCsl 


IOnC^ 


era' 


rvm-fr xj 


^ ^ ^ 04 CM 


O' O J 


c(\¥ r 


^3 >4T 


Si xj 


^ li 


1/ . Lf U ¥ U 


¥ 












¥ 


¥ 




* 








¥ 


O 




CJ 








a 



113 



r ajcrc>r-4fMm<j a a Or-t<\jro<^u or- cccrOf-^f\iro>t u or cccrOf-irvjfri^ inor- oocrOr-^fNirr 

ooc>f-^r-^r-:«-H,-Hf-H,-^r-^,-4f-4r\jcvj(\/rocvjroouc\jf>jrsjr.fnmfOfrif^rrjrofnrr!>r>J < ^ ^ -4 u^u^u^u^ir 

r r- r r* r r- r r- r- r-r r-^^f^r-r-^r-^ r-r--r- r- r- r- r- r-r r-h-r-h r-r-r-r r- r- r r-r-r r- 

fv c^lo^•c^lf^Jfv•rg^^)roc\lrof^^^^:f^Jf^Jf^lf^i^vif^Jf^♦^sJf^Jf^Jf^iror^Jc\^f^J(^:rgf^J^vJrsl(^J^o^^JCvc^^co<\jc^lc^jCMf^lC^4f^Jf^Jc^J 
UU.lLUUU ULlLlUUUULlLLLLUaaUUU.LLaU uuuu uuauuuuuu u uuuuuuuau 



•L 

•W 

* 

it 

V/ 

■yi- 

it 

it 

a 

cn 

U: 

> 

>— I 

K 

< 



U 

UJ 

o 

o 

2: 



ai 

if 

> 

-j 

u 

CL 

U'‘ 



IT 



LT* 

c 

o 



IP 

0 

1 

o 

o 

• 



IP 

rsj 



VP 

ro 



fp 



CO 





f- 








0 


1 




CD 




< 






a 


CO 






<I 




:;:? 






in IP 


2 


X 


oa. 


0 




•U 






r-a^ 


Z) 


Q 


CM«-H 


♦ 




< 






> *^3 


CO 


w 




'4- 




> 




CLO. 


OOtP 




CO 


•* •> 


1 




Ui 


tf. 


2 2 






cc ^ 


mo 


0 






0 


aa 




X 


<X -J 


mm 


0 




•1^ 






-J •‘X 




CJX 


•* #. 


• > 




it 


0 >xr»<\j 


•fr * 


+ Lr.o 


it 


z Ha 


om 


fM ■}<• 




♦ 








-> 


^ 2 OCL 




1 > 




* 


- _J>>-X> 


x>- 


0 ^ 


+ 


^ f-Hr— ( w- C^ 2 


r-H*— < 


— -r 






l-l w 1 + ♦ ♦ 


it it 


^ -- 1 


> 


0 II + a tt a 




>xx 


O' 


-fr 


+ II .-iu-1-i-Q.a 




xa> 






Cl-^^ 




m 


♦ 


^ >c*- + LJ X >• at. a. 


au 






V* 1 tt 


2H H 


UiXX 


t— ♦ 


a 


OUOOOw ♦ u u 


1 f\J(V 


II > ♦ xz)?~anr >-*>C3.-* ii 


lU 


OCD II II 


0 II x>* 0 


• • • • •a— 0 1 iixxKK 


x>r-r- 




1 


Q<io-> II w_ja 


H^uu a ct 


II <3<C 



K M M (I oooOf-^a:Dr-2:ci it it ^ 

•-^oca « r- ti It tt II It II i f\jfvj^> It ti II 

OUaCL^O>^>- I-KIIUO tl II i-hX 

* 



.ZZZJ 



IT. 



LTU) 

itUsO 

it 

it 

LJ 



m 

o 



X>C0COw 4 
II It ZTw It X 
l-I^U X II 
X >COh-»OX 

o in 

p- 



I +c»-^xau w2: xi-Hiii-HCLa 

II II 2:ouiu2^ It w.w.x>i-.w.cLX 

> ti acH uocH It xaa it it aa-jr> 
o>-»-iC5 



(DO 

r o 



U', 

o 



OU^ 0<\jix\ 
f-H ^ fNjrvjfvi 



114 



u vor-auo.~^<\jn>ju a cjr-iCNir -4 ir 

uM/ *i/ i* r^r* 

r-rr^r-r r r^r r r^r-r- r-r r-r r-r^ 

cvfMCvjfv rv.<\ rvr\toirvrvjr\tf\ rvj<Njf\jrgCMC\:f\fr\jCM 
UUUUUUU.UUUUUUUUUIXULLUL1U 
HKI- KKKl-Kh HKKI— 



I 



U. 

o 

o 

-» 

o 

w* 



OO 



UL 









c 








o 






W 








u 


<l 






o 


X 








o o 








—1 nC 






C" 


<3' *—1 






X 








a 


m 














<i 


•f r- 






•f 


o m 








w. o ST 






fNJirTLL 








-c XX p-i o 






ILfMO II II 


CJ •— 


o 


L‘ 


II II OUN 


-r 1 


• 


O-J— *Cs>f\J O •IT'O 


f— < pM 


ox 


II 




1 X 


It II 




O-f + II Z 


(vcox.2oo<i.-«NJu _;•— (xfv.::;:) cc:-~ct 


1 K II 


• t/ X 


1 




2: on xo> Q zru-^ 


U-Ca II --- 1 5i w»-CL 


II C II XX 


II X-I II OOC OOCMU. II I^ou^u 'Zi 


^ C X Z 2 <J 2T <_) O <_) Cf Q- 2" ^ •■ cr •-' a U^' 



O 

m 



U' 
r I 



ou^oio 
<r <» u tr 



o 


UAO 




nOO 




f— <fV’ 



115 



BIBLIOGRAPHY 



1. Thaler, G, J., Siljak, D. D. , and Dorf, R. C., 

Algebraic Methods for Dynamic Systems , University 
of Santa Clara, California, November I966. 

2. Hollister, F. H. , Network Analysis and Design by 

Parameter Plane Techniques , Ph.D. Thesis, U.S. Naval 
Postgraduate School, Monterey, California, 1965. 

3. Thaler, G. J. and Thompson, A. G., Parameter Plane 

Methods for the Study of the Frequency Response of 
Linear Time-Invariant Systems , paper presented at 
the 1967 Applied Mechanics Conference, Adelaide, 
Australia, June 196?. 

4. Glavls, G. 0., Frequency Response In the Parameter 

Plane , MS Thesis, U.S. Naval Postgraduate School, 
Monterey, California, I968. 

5. Staples, R. D. , Frequency Response Analysis of Non - 

linear Dynamic Systems , MS Thesis, U.S. Naval 
Postgraduate School, Monterey, California, I968. 



INITIAL DISTRIBUTION LIST 



No. Copies 
20 



1. Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 2231^ 



2. Library, Code 0212 2 

Naval Postgraduate School 

Monterey, California 939^0 

3. Naval Ship Systems Command (Code 2052) 1 

Department of the Navy 

Washington, D. C. 20360 

Prof. Geo. J. Thaler 10 

Department of Electrical Engineering 
Naval Postgraduate School 
Monterey, California 939^0 

5. Prof. S. R. Parker 1 

Department of Electrical Engineering 

Naval Pobtgraduate School 
Monterey, California 939^0 

6. LCDR R. A. Nichols, USN 1 

USN Underwater Sound Laboratory 



Port Trumbull, New London, Connecticut 06320 



117 



UNCLASSIFIED 

Security Classification 



DOCUMENT CONTROL DATA - R & D 



fSerurify' classification of title, body of abstract and indcxinti annotation must be entered when the overall report i» c laMnified) 



\ ORIGINATING A,CT\^\ T y (Corporate author) 



2a. REPORT SECURITY C L A S 51 P I C A T I ON 



Naval Postgraduate School 
Monterey, California 939^0 



Unclassified 

2b. CROUP 



3 REPORT TITLE 



PARAMETER PLANE DESIGN OF ELECTRICAL 
FILTERS CONTAINING NONLINEAR ELEMENTS 

4 DESCRIPTIVE NOTES (Type of report and^^inclusive dates) 

Electrical Engineer's Thesis; December 1969 

5 Au THO R(S> ('Firsr name, rriiddle initial, last name) 

Roger Alan Nichols, Lieutenant Commander, United States Navy 



C REPOR T D A TE 

December 1969 

Ba. CONTRACT OR GRANT NO. 



7a. TOTAL NO. OF PAGES 176. 

1]J |_ 

9a. ORIGINATOR'S REPORT NUMBER 



NO 



(S) 



OF REFS 

5 



6. PROJEC T NO. 



96. 



OTHER REPORT NO(S) (Any of6cr numbers that may be aa st ^ed 
this repor(^ 



d. 



10, DISTRIBUTION STATEMENT 



This document has been approved for public release and sale; 
its distribution is unlimited. 



11. SUPPLEMEN T ARY NOTES 



12. SPONSORING MILITARY ACTIVITY 



Naval Postgraduate School 
Monterey, California 939^0 

13 ABSTRACT ^ p ar amo t o r~ p 1 an o method, based upon a quasi-frozen 

assumption, for the design of electrical filters containing 
slowly varying nonlinear elem.ents is presented together 
with examples of its application to the lowpass, highpass, 
bandpass and band-reject filter types. The nonlinear mag- 
nitude and frequency characteristics are presented in the 
form of parameter-plane graphs which allow value judgements 
to be made concerning the choice of nonlinear elements used 
in the filter design. The computer programs used to generate 
the graphs are also Included. 



DD "o'.“..1473 (PAGE 1) 

S/N 0101 -807-681 I 



119 



UNCLASSIFIED 

Security Classification . . 

A-31400 




UNCLASSIFIED 

Security Classification 



1 4 

KEY WORDS 


LINK A . 


LINK B 


LINK C 


ROLE 


WT 


ROLE 


W T 


ROLE 


W T 


ELECTRICAL FILTER 
NONLINEAR ELEMENTS 
PARAMETER PLANE 
COMPUTER GRAPHICS 















DD .1473 UNCLASSIFIED 



S/N 0101 - 807-6821 



120 



Security Classification 



